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Abstract 

We present new CCD photometry of the distant old open star cluster Berkeley 32 in Johnson V 
and Cousins I passbands. A total of ~ 3200 stars have been observed in a field of about 13 x 13 . 
The colour-magnitude diagram (CMD) in V, (V — I) has been generated down to V = 22 mag. 
A broad but well defined main sequence is clearly visible. Some blue stragglers, a well developed 
subgiant branch and a Red Clump are also seen. By fitting isochrones to this CMD as well as to 
I—} | other CMDs available in the literature, and using the Red Clump location, the reddening, distance 

and age of the star cluster have been determined. The cluster has a distance of ~ 3.3 kpc, its radius 
is about 2.4 pc; the reddening E(B-V) is 0.08 mag and the age is ~ 6.3 Gyr. By comparison with 
7—1 . theoretical isochrones, a metallicity of [Fe/H] « —0.2 dex has been estimated. 

Theoretical isochrones have been used to convert the observed cluster luminosity function into 
| a mass function in the mass range ~ 0.6—1.1 Mq. We find a much flatter mass function than what 

has been found for young clusters. If the mass function is a power law dN ~ m a dm, then we get 
a = —0.5 ± 0.3. This may be seen as a signature of the highly evolved dynamical state of the 
I . cluster. 
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to ! 1 Introduction 

i Berkeley 32 (C0655 + 065 - OCL 522, 1=207.°95, b=4.°4), also known as Biurakan 8, is 
^ ! a small galactic (open) star cluster of angular diameter ~ 6'. It is located in the Galactic 
5^ ! anticentre direction and has been classified as Trumpler class II2m (Lynga 1987). This object 
was discovered by Iskudarjan (1960) and catalogued by Setteducati & Weaver (1960). On 
the sky survey maps, the cluster appears to be rich and likely of old age (cf. King 1964). 
The first photometric study of the cluster was carried out by Kaluzny & Mazur (1991) in the 
UBV and Washington systems by CCD imaging of an area ~ 6. 6 x 6. 6. They presented the 
colour-magnitude diagram (CMD) and discussed its morphology, derived cluster reddening and 
distance as E(B — V) = 0.16 mag and 3.1±0.2 kpc respectively, and estimated the metallicity as 
[Fe/H] = —0.37 ±0.05. They derived an age of ~ 6 Gyr. Using morphological age parameters, 
Janes & Phelps (1994) estimated an age of 7.2 Gyr for Berkeley 32. 

Scott et al. (1995) have determined radial velocities of 7 cluster members and found a mean 
heliocentric radial velocity of +101±3 km/sec, which among their cluster sample deviates most 
from a pure Galactic rotation after Berkeley 17. 

Mass function studies of open star clusters indicated that the slopes of the mass functions 
of older (age > 1 Gyr) clusters differ significantly from each other and also are not uniform 
over the entire mass range (cf. Sagar & Griffiths 1998b). However, the number of old objects 
studied so far is small. As such, a study of the old cluster Berkeley 32 can contribute interesting 



Our aim is to re- analyse Berkeley 32 with deeper photometry than Kaluzny & Mazur (1991) 
had at their disposal and study its mass function which is lacking. The observations, data 
reductions and comparison with earlier photometry are given in the next two sections. The 
cluster radius, other photometric results and mass function of the cluster are described in the 
subsequent sections of the paper. 



2 Observations and Reductions 

The CCD observations have been obtained on 1999 March 22 with the 3.5m telescope at Calar 
Alto Observatory, Spain, run by the Max-Planck Institute for astronomy, Heidelberg. The focal 
reducer MOSCA attached at the RC Cassegrain focus provided an effective f/2.7 focal ratio 
flhttp: / / www. caha.es / caha / instruments / mosca / manual. html ) . The observations have 



been carried out in the Johnson V and Cousins I filters using the SITE 18b CCD chip. Each 
pixel (24// square) corresponded to 0."53 x 0. '53 on the sky. The non-vignetted area of the 
CCD was 1500 x 1500 pixel 2 providing a field of about 13. 25 x 13. 25. The read-out noise was 
5.4 electrons per pixel and the ratio electrons-per-ADU was ~ 2.7. Figure 1 shows the finding 
chart for the imaged cluster region and Table 1 lists the log of the observations. For calibration 
purpose, we observed the standard star field SA 98 (Landolt 1992). 

The present observations were carried out as a back-up programme. We therefore could 
not observe standard stars during the whole night. Instead, we observed the standard field 
SA 98 three times: before, between, and after the observations of Berkeley 32 (Table 1). The 
strategy was to observe Berkeley 32 at a similar air-mass as the standard fields, so that one 
can calibrate photometric data of the cluster region with an accuracy of a few percent without 
a precise determination of the atmospheric extinction coefficients. The basic processing of the 
data frames was done in the standard manner using the MIDAS data reduction package. The 
uniformity of flat fields is better than one percent in both filters. 

Instrumental magnitudes have been measured using the DAOPHOT software (Stetson 1987, 
1992) under MIDAS. The image parameters and errors provided by DAOPHOT were used to 
reject poor measurements. About 10% of the stars were rejected in this process. In those cases 
where brighter stars are saturated on deep exposure frames, their magnitudes have been taken 
only from the short exposure frames. Most of the stars brighter than V ~ 12 mag could not 
be measured because they were saturated even on the shortest exposure frames. The CCD 
instrumental magnitudes have been calibrated using the observations of the SA 98 field and 
the following relations 

V - v = aO + al * (V - I); (V - I) - (v - i) = bO + bl * (V - I) 

where capital letters denote standard magnitudes and colours, and lower case letters denote 
instrumental values. The values refer to exposure time of 1 second. These equations along with 
the site mean atmospheric extinction values of 0.15±0.04 and 0.09±0.02 mag per unit air-mass 
in V and (V — /) respectively were used in determining the colour equations for the system 
as well as the zero-points. The effects of uncertainties in atmospheric extinction values are 
maximum on zero-points but least on the colour coefficients. We therefore averaged the colour 
coefficients from the individual standard observations. With these values fixed we calculated 



Table 1. Log of CCD observations of the cluster Berkeley 32 and the calibration region SA 98 (Landolt 
1992). The data have been obtained on March 22, 1999. 



Object 


Time 


Filter 


Exposure time 


Air-mass 




(UT) 




(seconds) 




SA 98 


19:23:45 


V 


8 


1.27 


SA 98 


19:29:00 


I 


8 


1.27 


Be 32 


19:40:04 


v 


10 


1.17 


Be 32 


19:44:57 


V 


60 


1.18 


Be 32 


19:49:05 


I 


8 


1.18 


Be 32 


19:51:47 


I 


40 


1.18 


Be 32 


19:54:36 


I 


3 


1.19 


SA 98 


20:00:04 


V 


10 


1.30 


SA 98 


20:02:37 


I 


8 


1.30 


Be 32 


20:05:29 


V 


10 


1.20 


Be 32 


20:07:44 


V 


100 


1.20 


Be 32 


20:11:37 


I 


8 


1.21 


Be 32 


20:14:04 


I 


80 


1.21 


SA 98 


20:43:19 


V 


10 


1.39 


SA 98 


20:45:37 


I 


8 


1.40 



Table 2. Colour coefficients and zero-points are for 1 second exposure time of the standard stars. 



Air-mass 


a0±cr 


al±cr 


b0±C7 


bl±cr 


1.27 


-0.513 ±0.01 


0.07 ±0.01 


0.667 ±0.02 


0.03 ± 0.003 


1.30 


-0.523 ±0.02 


0.07 ±0.01 


0.664 ± 0.02 


0.03 ± 0.003 


1.39 


-0.573 ±0.02 


0.07 ±0.01 


0.643 ± 0.02 


0.03 ± 0.003 



Table 3. Internal photometric errors as a function of brightness, a is the standard deviation (a) per 
observations in magnitude. 



Mag 


range 


ay 


07 


12.0 


- 14.0 


0.005 


0.010 


14.0 


- 16.0 


0.005 


0.010 


16.0 


- 17.0 


0.005 


0.010 


17.0 


- 18.0 


0.006 


0.013 


18.0 


- 19.0 


0.009 


0.024 


19.0 


- 20.0 


0.017 


0.051 


20.0 


- 21.0 


0.041 


0.117 


21.0 


- 22.0 


0.106 




22.0 


- 23.0 


0.221 





zero-points are uncertain by ~ 0.02 mag in V and (V — I). The internal errors as a function 
of magnitude for each filter are given in Table 3. The errors become large (> 0.15 mag) for 
stars fainter than V = 22 and I = 21 mag. The X and Y pixel coordinates as well as the 
V and (V — I) magnitudes and DAOPHOT errors of the stars observed in Berkeley 32 are 
listed in Table 4. Stars observed by Kaluzny & Mazur (1991) have been identified in the last 
column. Table 4 is available only in electronic form at the open cluster database Web site at 
[http:/ /obswww.unige.ch/webda/| . It can also be obtained from the authors. 



3 Comparison with previous photometry 

In this section, we compare the present CCD photometry with the only previous CCD photo- 
metric observations of the cluster by Kaluzny & Mazur (1991) in the only common passband V. 
The transformation equations relating their (Xk m , Yk m ) coordinate system to ours (X pres , Y pres ) 
were found to be 

X pres = 1178.252 - 0.024X fcm - 1.512F fcm ; Y pres = 259.383 + 1.511X fcm - 0.024F fcm 

There are 835 stars measured by Kaluzny & Mazur (1991) whose positions coincide within 1 
pixel with the stars positions measured by us. The differences (AV) between the present data 
and data obtained by them are plotted in Fig. 2, while the statistical results are given in Table 
5. These show that except for a few outliers, which appear to be mostly stars that were treated 
as single in one study and as double (due to blending) in the other, the distribution of the 
photometric differences seems fairly random with almost no zero-point offset. As expected, the 
scatter increases with decreasing brightness and becomes more than ~ 0.1 mag at fainter levels. 
Considering the uncertainties present in our and Kaluzny & Mazur's (1991) measurements, we 
conclude that they are in very good agreement. 



Table 5. Statistical results of the photometric comparison with data from Kaluzny & Mazur (1991). 
The difference (A) is in the sense present minus comparison data. The mean and standard deviation 
(cr) are based on N stars. A few points discrepant by more than 3.5 a have been excluded from the 
analysis. 



V range 


AV 




(V — I) range 


AV 




(mag) 


Meanicj 


N 


(mag) 


Meanicr 


N 


12-14 


0.006±0.03 


17 


-0.1-0.65 


-0.021±0.16 


41 


14-16 


-0.012±0.04 


67 


0.65-0.8 


-0.013±0.12 


197 


16-17 


-0.004±0.04 


129 


0.8 -1.0 


-0.020±0.12 


242 


17-18 


-0.015±0.08 


144 


1.0 -1.5 


-0.026±0.11 


286 


18-19 


-0.029±0.09 


159 


1.5 -3.2 


-0.011±0.11 


69 


19-20 


-0.015±0.14 


157 








20-21 


-0.023±0.16 


162 









4 Radius of the cluster 

We used radial stellar density profile for the determination of cluster radius. Such determina- 
tions can provide ambiguous results as it depends on the limiting magnitude of the star counts. 
The fainter the stars are, the larger becomes the cluster radius, if mass segregation due to 
two-body relaxation is present. Given these caveats, it is not our aim to derive a dynamically 
relevant radius, but to determine the region where the cluster population dominates over field 
stars so that it can be used for investigations of the cluster properties. 

We derived the position of the cluster centre by iteratively calculating the average X and Y 
positions of stars within 150 pixels from an eye estimated centre, until it converged to a constant 
value. The (X,Y) pixel coordinates of the cluster centre are (700, 665) with an accuracy of few 
pixels. The radial stellar density profile determined up to ~ 6' from the cluster centre using 
stars brighter than V=18 mag is plotted in Fig. 3. The radius at which the star density flattens 
is considered as cluster radius which is ~ 2.' 7. This agrees fairly well with the value of 3' given 
by Lynga (1987). We fit the following form of a King (1962) profile to the observed stellar 
density distribution 

/(r)ocC-(^ = ) 2 , 

where C is the central stellar density, r c and r t are the core and tidal radius respectively. A 
least square fitting of the profile to the observed points yielded C = 33.9 ± 8 star/ arcmin 2 , r c = 
1.0 ± 0.38, r t = 23 ± 50 arcmin. 



5 Colour-magnitude diagrams (CMDs) 

5.1 The V, (V - I) CMD 

We plot the V, (V — I) CMD for all (~ 3200) measured stars in the region of Berkeley 32 in 
Fig. 4 (A). The CMD reaches down to V = 22 mag. The cluster main-sequence (MS) contam- 



region is of composite nature, the cluster population appears to be dominating. The only way 
to sharpen morphological features of the cluster sequence in the CMD is to select stars with 
small radial distances by compromising between a decreasing number of cluster stars and an 
increasing field population. The Fig. 4 (B) shows our best result. Here we have selected only 
stars with a radial distance up to ~ 2. '7. The features of a very old open star cluster namely 
the distinct turn-off region and the subgiant branch are now very clearly visible. The giant 
branch (GB) is very sparsely populated and not well defined. Moreover, a group of stars can be 
seen which are brighter and bluer than the MS turn-off point suggesting that some of them are 
blue stragglers (BSs). Such objects have been found in most intermediate and old age open star 
clusters (see Kaluzny 1994; Phelps et al. 1994; Sagar & Griffiths 1999a). Many of them have 
been identified as close binary systems. The (X,Y) pixel coordinates, radius, magnitudes and 
colour of the stars located in the GB, red GB and BS regions of the CMDs are given in Table 
6. The cluster membership of these stars is also indicated in the table. A star is considered 
as probable cluster member if it lies within ±0.05 mag in colour and ±0.1 mag in brightness 
with respect to the isochrone of the cluster age at least in two of the V, (U — V); V, (B — V) 
and V, (V — I) diagrams. In addition, brightening due to unresolved/optical binary stars has 
also been considered. 

Fig. 4 (C) shows the V, (V — I) diagram of stars with radial distances more than ~ 4. '4 from 
the cluster centre. Overplotted are the fiducial points of the cluster sequence. There are a few 
red giants, which perhaps still belong to the cluster. A considerable part of the main sequence 
population has a turn-off similar to the cluster, but the bulk of the main sequence stars are 
clearly shifted towards the red, indicating higher reddening, and thus a background population. 
However, the interesting question whether there are evaporated stars surrounding the cluster 
can not be answered on the basis of the present data. For this, kinematic informations like 
proper motions and radial velocities of these stars are required. 

5.2 The cluster age from the "Red Clump" 

It is well known that for intermediate and old open star clusters, the location of the Red Clump 
(RC) (the more massive analog of the horizontal branch in globular clusters) relative to the MS 
turn-off point is correlated with age (cf. Kaluzny 1994; Phelps et al. 1994; Carraro & Chiosi 
1994; Pandey et al. 1997 and references therein). The two morphological parameters generally 
used for estimating cluster ages are the differences in magnitudes ( AV) and colours {A(B — V) 
or A(V — I)) between the RG branch at the level of the clump and the MS turn-off point, with 
the advantage that no prior knowledge of cluster distance, reddening and accurate metallicity 
is required. 

Following Kaluzny (1994), we find AV = 2.7 ± 0.05, A(V - I) = 0.45 ± 0.03 in the case 
of Berkeley 32. Using the relation given by Carraro & Chiosi (1994), we derive log (age) = 
9.8 ±0.1 for the cluster. A slightly modified version of AV has been introduced by Janes 
& Phelps (1994) who used the luminosity difference between the RC and the inflection point 
between the turn-off region and the subgiant branch. We confirm their value of 2.4 for Berkeley 
32. From their Fig.l one reads off the logarithm of age as 9.8-9.9, yielding an age of 6.3 - 8 
Gyr. According to their Table 1, there are only a few clusters older than Berkeley 32, like NGC 
6791, Berkeley 54, AM 2 and Cr 261. It is thus clear that Berkeley 32 belongs to the group of 



very old open clusters in our Galaxy. 

5.3 Determination of the cluster parameters using theoretical isochrones 

We have determined the colour excess, the distance modulus, and also the age of the cluster 
by fitting theoretical stellar evolutionary isochrones from the set of Bertelli et al. (1994) to our 
V, (V — I) diagram. These isochrones are derived from stellar models computed with updated 
radiative opacities and include the effects of mass loss and convective core overshooting. The 
models trace the evolution from the zero-age main-sequence (ZAMS) to the central carbon 
ignition for massive stars and to the beginning of the thermally pulsing regime of the asymptotic 
giant branch phase for low and intermediate mass stars. 

As most of the factors responsible for the colour spread in the MS will redden the stars 
(differential reddening, binaries, rotation, star spots), we have used the blue envelope of the 
MS in the CM diagram for the estimation of the cluster parameters. We fit the isochrones by 
eye taking into account the observational error. It turns out that the isochrone with log(age) 
= 9.8, X = 0.7, Y = 0.28 and Z = 0.008 fits best to the cluster locus, including the RC, and 
thus is in good agreement with the age estimated from the morphological parameters of the 
cluster CMD. In order to also demonstrate upper limits of the effects of binaries in the CMD, 
the log (age) = 9.8 isochrone for the single stars has been brightened by 0.75 mag leaving the 
colour unchanged. A maximum reddening of E(V — I) = 0.11 (or E(B — V) — 0.08) mag 
can be applied to place the isochrone correctly on the cluster sequence observed in Fig. 4(B). 
Some stars above the turn-off point lie on the isochrones of binaries indicating the possibility of 
being indeed binary members of the cluster. The lower giant branch in the V, (V — I) diagram 
appears marginally too blue, indicating that the cluster may have a slightly higher metallicity. 
However, a Z = 0.02 isochrone is definitely too metal-rich. Moreover, a solar metallicity would 
decrease the cluster reddening even further, while Kaluzny & Mazur (1991) quote a reddening 
of E(B — V) — 0.16 mag. However, such a high reddening is supported only by their V, (U — V) 
diagram (see section 5.4). 

It can also be seen that the theoretical location of the RC fits rather well with the observed 
one for Berkeley 32 unlike in some other old open star clusters (see Sagar & Griffiths 1999a). 
For example, it is too faint for NGC 6603 and too bright for NGC 7044. 

The value of the apparent distance modulus derived from Fig. 4(B) is 12.8 mag. Here we 
adopt the reddening law of Rieke & Lebofsky (1985), who give Ay/E(V — I) = 1.94. With 
the extinction Ay = 0.21 mag, we get for Berkeley 32, a true distance modulus of 12.6 with 
an uncertainty of ~ 0.15 mag which includes errors in the photometric calibration, isochrone 
fitting and the reddening determination. 

5.4 Isochrone fitting to the UBV data of Kaluzny & Mazur (1991) 

Fig. 5 shows the V, (U — V) and V, (B — V) diagrams generated from the Kaluzny & Mazur 
(1991) photometric data. Overplotted are the isochrones of Bertelli et al. (1994), having the 
same ages, metallicity and helium abundance as we used for our V, (V — I) diagram in Fig 
4(B). In order to fit the isochrone to the cluster sequence, we had to employ a reddening of 
E(B -V) = 0.08 mag and E(U - V) = 0.22 mag in the V, {B - V) and V, (U - V) diagrams 
respectively. While for the reddening law of Rieke & Lebofsky (1985), the E(B — V) value 



Table 7. The age, metallicity ([Fe/H]), reddening (E(B - V) and E(V - I)), true distance 
modulus((m — M)o), distance, galacto-centric distance (Rgc) (adopting galacto-centric distance of 
the Sun as 8 kpc) and z-distance for Berkeley 32. 

Age 6.3 Gyr 

[Fe/H] -0.2 dex 
E(V — I);E(B — V) 0.11;0.08 

(m-M)o 12.60 ±0.1 

Distance 3.3 kpc 

R GC 11.0 kpc 

z 250 pc 



agrees well with our E(V-I) value (E(V - I)/E(B - V) = 1.6), the E(U - V) is too large 
(E(U - V)/E(B - V) = 1.6A). On the other hand, it is too small for E(B - V) = 0.16 mag, 
given by Kaluzny & Mazur (1991) and is thus not compatible with the reddening values derived 
from the V, (V—I) and V, (B—V) diagrams. This may suggest that the [/-photometry is perhaps 
in error and we adopt E(B — V) = 0.08 mag as the value for the cluster reddening. 

5.5 The cluster distance from the Red Clump 

For a star cluster as old as Berkeley 32, an attractive method to determine its distance is using 
the location of the RC of intermediate-age helium core burning stars as a standard candle (e.g. 
Paczynski & Stanek 1998). The absolute I-magnitude of RC stars in the solar neighborhood has 
been calibrated by Hipparcos parallaxes, resulting in Mf = —0.23 ±0.03. Cole (1998) discusses 
the age and metallicity dependence of the RC-brightness and notes that for poulations older 
than 4-5 Gyr, the M° is independent of stellar mass, but still shows a metallicity dependence 
of the order 5Mf = (0.21 ± 0.07) [Fe/if], where 5 Mi is the brightness difference between the 
RC in the solar neighborhood and the population under consideration. As the value is small 
for Berkeley 32, we neglect this correction here. 

The RC in Berkeley 32 has V = 13.67 ± 0.03 and (V - I) = 1.16 ± 0.03, where the error is 
the uncertainty in the definition of the RC in the CMD. This yields (m — Mj) = 12.74 ± 0.08 
as the apparent distance modulus, if we include the photometric calibration uncertainty in the 
error. The extinction in I is determined by Aj/E(V — I) — 0.93 which is 0.10 for a value of 
E(V — I) — 0.11. If we assign an additional error of 0.05 to the extinction correction, we have 
(m — M) = 12.64 ±0.1 as the value for true distance modulus. This agrees well with the value 
obtained using isochrone fitting. In the following, we therefore adopt 12.6 ±0.1 as the value for 
the distance modulus of the cluster. The present distance determination of 3.3±0.2 kpc agrees 
very well with the value of 3.1 kpc given by Kaluzny & Mazur (1991). The cluster parameters 
derived by us are listed in Table 7. 

5.6 Location of Berkeley 32 in the Galaxy 

The cluster Berkeley 32 occupies an important position for understanding the variation of 
metallicity in the Galactic disk, as the issue of the existence of a metallicity gradient is not 
yet settled. According to Friel (1995), the metallicities of open clusters indicate a gradient of 



—0.09 dex/kpc. On the other hand, Twarog et al. (1997) argue that the open cluster system 
can be divided in 2 radial groups, with a very flat or even vanishing gradient in each group. 
Their mean metallicities differ by 0.3 dex and there is a discontinuity at a radial distance of 10 
kpc. As the galactocentric distance of Berkeley 32 puts it just near this discontinuity, Berkeley 
32 could help to decide between these two metallicity patterns in the Galactic disk. However, a 
more accurate determination of the metallicity than we are able to do, is required. Also, more 
clusters/objects either in the vicinity of Berkeley 32 or at similar galacto-centric distances need 
to be observed before the the metallicity pattern can be unambiguously determined. 

6 Mass function 

The study of the mass function (MF) of Berkeley 32 is based on a pair of deep V and / CCD 
frames only. This is done for evaluating the data completeness accurately. Guided by the radial 
stellar density profile in Fig. 3, we selected stars located within a circle of 165 arcsec radius 
(surface area 23.76 square arcmin) around the cluster centre for the MF study. With the aim 
of detecting possible radial MF variations, the data completeness has first been evaluated in an 
inner and outer region separately, but it turned out to be the same within the errors. Besides, 
the small number statistics prevented us from a detailed study. We therefore determined the 
MF for the entire region without any subdivision. 

To suppress the field star contamination as far as possible, we selected cluster main sequence 
stars by using the following boundaries in the V, (V — I) diagram: 

V > 8.91 + 9.45 -(V-I) + 1.8 • (V - I) 2 

and 

V < 11.52 + 12.64 .(y-I) + 3.86 • (V - I) 2 

The field star contamination has been determined using the remaining chip area outside a radius 
of 265 arcsec (surface area 123.1 square arcmin) from the cluster centre. Table 8 lists the field 
and cluster counts derived in this way along with their completeness factors. The completeness 
factors have been determined by using artificial stars along the clusters main sequence and 
recovering them in the CMD, not in one filter, to be as realistic as possible. Table 8 also lists 
the numbers which are relevant for the MF determination. The transformation from apparent 
to absolute visual magnitude (My) has been done using the cluster parameters given in Table 
7. The isochrone log(age) = 9.8 and z = 0.008 provides the following parametrization of mass 
(m) and My: 

m = 1.665 - 0.186 • My + 0.00698 • M v . 

The values of the normalised counts (N) are in stars/arcmin 2 . They are corrected for com- 
pleteness and field star contamination, and divided by the mass interval of the magnitude bin. 
The errors of the normalised counts result from error propagation. This may be incorrect from 
a puristic viewpoint, as they are no longer small compared to the counts. However, we do not 
use them any further beyond a qualitative demonstration that they are large. 

A common description of the stellar mass spectrum is a power law dN oc m a dm, where dN 

f +ko A/TT7 TT 



Table 8. The V-magnitude of the bin center, the raw counts for cluster (Nc) and field regions (Np), 
and the corresponding completeness factors fc and fp are listed. Absolute My and stellar mass (m) 
of the bin center, the mass interval (Am) corresponding to the magnitude bin, the normalised counts 
(TV) and their errrors (SN) are also listed. 



V 


N c 


N F 


fc 


fF 


My 


m 


Am 


N 


8N 


16.25 


41 


53 


1.00 


1.00 


3.40 


1.113 


0.0693 


18.70 


4.2 


16.75 


44 


67 


1.00 


1.00 


3.90 


1.046 


0.0658 


19.88 


4.6 


17.25 


47 


86 


1.00 


1.00 


4.40 


0.982 


0.0623 


20.54 


5.0 


17.75 


43 


107 


1.00 


1.00 


4.90 


0.921 


0.0588 


16.00 


5.1 


18.25 


50 


135 


0.95 


1.00 


5.40 


0.864 


0.0553 


20.22 


6.3 


18.75 


37 


121 


0.85 


1.00 


5.90 


0.811 


0.0518 


16.39 


6.3 


19.25 


35 


127 


0.80 


0.93 


6.40 


0.761 


0.0483 


15.15 


7.4 


19.75 


31 


158 


0.81 


0.89 


6.90 


0.714 


0.0449 


3.76 


7.4 


20.25 


37 


154 


0.80 


0.90 


7.40 


0.671 


0.0413 


13.46 


8.8 


20.75 


31 


142 


0.73 


0.86 


7.90 


0.631 


0.0379 


11.73 


8.9 
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a uniform exponent is at best realised within limited mass intervals. The universality of the 
slope of the initial mass function (IMF) is still a matter of discussion (for a recent review see 
Scalo 1998), but studies of a large number of young clusters in the Milky Way and the Large 
Magellanic Clouds do not speak evidently against an universal IMF at least above 1 M (e.g. 
Sagar et al. 1986; Sagar & Richtler 1991; Janes & Phelps 1994; Fischer et al. 1998; Sagar 
2000), with a around —2.3. Below 1 M , the data for young open clusters are sparse and any 
secure statement on a possible universal IMF is not yet possible. 

However, in the case of a very old cluster like Berkely 32, we might anyway not expect the 
IMF to be still realized. As a cluster evolves dynamically, low mass stars evaporate out of the 
cluster potential faster than high mass stars. In a cluster much older than its relaxation time, 
the dynamical effect therefore can change an originally rising IMF into a flat or even declining 
MF. 

Fig. 6 shows the MF for Berkeley 32. The logarithm of mass is plotted against the logarithm 
of the normalised counts. Note that the binning in mass is also logarithmic. Although the errors 
are too large for any deeper analysis, it is apparent that the MF is much natter than of most 
young clusters. A fit to a power-law indeed returns the value a = —0.5 ±0.3 while we expect in 
this mass domain an exponent around —2 for young clusters (Richtler 1994). Here it must be 
remarked that we assumed single stars for the applied mass-luminosity relation of the isochrone. 
Sagar & Richtler (1991) discussed how the presence of binaries flattens a "true" IMF. But even 
if we find a large binary fraction in Berkeley 32, as the CMD suggests, their effect would by far 
not be sufficient to steepen the observed MF to the level as it is observed for young clusters. 

This behaviour, in agreement with theoretical expectations, has been found for other old 
open clusters as well. For example, Francic (1989), among his sample of 8 open clusters, found 
the old objects NGC 752 (2.5 Gyr) and M67 (5 Gyr) to show even declining mass functions 
with a > 0. Further well studied open clusters like NGC 6791 (Kaluzny & Rucinski 1995), 
NGC 188 (von Hippel & Sarajedini 1998), NGC 2243 (Bergbusch et al. 1991) also have flatter 



MFs. But one also can find old clusters with MF not distinguishable from a Salpeter mass 
function, e.g., mass spectrum of Berkeley 99 (age 3.2 Gyr) has a ~ —2.4 (Sagar & Griffiths 
1998b). This demonstrates that open clusters do have distinctly different dynamical histories, 
which may depend on their structure, total mass, location, orbit characteristics etc. 

7 Summary 

New V and / CCD photometry down to V = 22 mag is presented for about 3,200 stars in 
the region of the open cluster Berkeley 32. The present photometry serves as a data base for 
determining the cluster properties and to study the stellar mass function for the first time. 
The cluster's radial density profile is well represented by a King (1962) profile. By fitting 
of theoretical isochrones and using the location of the Red Clump, we confirm earlier results 
that it is indeed a very old open cluster (6.3 Gyr). Its metallicity is between Z = 0.008 and 
Z = 0.02, distance is 3.3 kpc and galacto-centric distance is 10.8 kpc. Clusters/objects with 
these characteristics can play a very valuable role to distinguish between the two models of the 
metallicity variation in the Galactic disk, advocated by Friel (1995) and Twarog et al. (1997) 
respectively. However, the case of Berkeley 32 is ambiguous. The parameters of Berkeley 32 are 
compatible with both a smooth Galactic metallicity gradient as well as with its membership of 
the cluster population of the inner domain of Twarog et al. (1997). 

We also investigated the mass spectrum of Berkeley 32 in the mass range 0.6-1.1 M Q . A 
power-law fit returns a = —0.5 ± 0.3 for the slope of the MF, which is much flatter than the 
slopes found in young open clusters. Berkeley 32 shares this behaviour with other old open 
clusters which indicates an evaporation of low-mass cluster stars. 
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Figure captions 

Fig. 1 Identification chart for the Berkeley 32 region. The (X,Y) coordinates are in CCD pixel 
units and one CCD pixel corresponds to 0."53 on the sky. North is up and East is to the 
left. Filled circles of different sizes represent the brightness of the stars. The smallest size 
denotes stars of V = 17 mag. 

Fig. 2 Comparison of the present V magnitude with those of Kaluzny & Mazur (1991). The 
differences (A) are in the sense of this study minus Kaluzny & Mazur. They are plotted 
against the present CCD photometry. 

Fig. 3 Plot of the radial density profile (•) for stars brighter than V = 18 mag in Berkeley 
32 region. The length of the bar represents errors resulting from sampling statistics. 
Overplotted (solid curve) is a King (1962) profile with parameters given in the text. The 
arrow denotes the radius where the surface density of cluster stars becomes becomes 
comparable with the field star density. 

Fig. 4 The V, (V — I) diagrams (A) for all stars observed by us, (B) for the cluster population 
(stars with radius < 2. '7) and (C) for the field population (stars with radius > 4/4) in the 
Berkeley 32 region are plotted. The composite nature of the stellar population is apparent 
in (A). In the cluster population, the Bertelli at al. (1994) isochrones for z = 0.008 and 
log (age) = 9.7 (short-dashed), 9.8 (continuous) and 9.9 (dot-short dashed) are shown. 
The isochrone of log (age) = 9.8 found to be best fitting to the observed cluster sequence 
with a reddening of E(V — I) = 0.11 (or E(B — V) — 0.08) and an apparent distance 
modulus of 12.8 mag. The dotted curve shows the extent that binaries of equal mass can 
brighten the isochrone of log (age) = 9.8. In the field population, the overplotted curve 
is displaying the cluster locus. There are red giants resembling cluster giant branch stars. 
A considerable part of the main sequence population has a turn-off similar to the cluster, 
but the bulk of the main sequence stars are clearly shifted towards the red, indicating 
higher reddening, and thus a background stellar population. 

Fig. 5 The V, (U — V) and V, (B — V) diagrams generated for the cluster population of Berkeley 
32 from the photometric data of Kaluzny & Mazur (1991). The Bertelli et al. (1994) 
isochrones of the same metallicity and ages, shifted by the same value of the apparent 
distance modulus as in Fig. 4(B) are shown. The best eye fits to the cluster sequence 
for the reddening values shown in the plot. The E(U — V) value is not compatible with 
E(B — V), indicating a problem with the U photometry. 

Fig. 6 The mass function of Berkeley 32 between 0.6 and 1.1 M Q as derived from the present 
data. Although the error bars are large, the mass function is clearly flatter than what has 
been found for young star clusters. 



Table 6. Spatial and UBVI photometric values of the candidate giant branch and blue straggler 
stars are listed along with the identification of Kaluzny & Mazur (1991) prefixed with KM. The 
(U — B) and (B — V) values are taken from Kaluzny & Mazur (1991). The probable photometric 
members have been identified as PM in the last column. 



(A) Stars redder than MS turn-off (candidate for giant branch) 
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(B) Stars bluer than the MS turn-off (candidate for blue straggler) 
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